Gapless Superfluid in SU(2n+ 1) Fermions Systems 
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We investigate the generalized Hubbard model of (2n + 1) Fermion species interacting via a 
symmetric contact attraction potential. We prove that the ground state of such system is a gapless 
superfluid, where a full Fermi surface coexists with a superfluid. Moreover, doing so we prove the 
existence of a free mode in a strongly interacting system, regardless of the potential strength. This 
proof holds at the mean field level. A Grassmannian Gaussian integration technique is used to deal 
with the problem. Our predictions may be relevant to future high-spin cold atoms experiments. 
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The smooth crossover from a Bardeen-Cooper- 
Schrieffer (BCS) superfluidity, valid for weakly interact- 
ing Fermions, to a Bose- Einstein Condensation (BEC), 
valid for (composite) Bosons, was suggested theoretically 
about 30 years ago jl|. This prediction was verified ex- 
perimentally in the last decade by several groups using 
trapped cold atoms. Two-component Fermionic atoms, 
like 40 K [2] or 6 Li were trapped and the Feshbach 
resonance was used to control the interaction between 
them. 

Recently, higher spin atoms were trapped, for which F, 
the hyper-fine spin of the atom, enables studying systems 
of (2F + 1) Fermion species. For example, an F = 1 
system was materialized with 6 Li atoms [4[, an F = 5/2 
with 173 Yb atoms and F = 9/2 with 87 Sr atoms @. 
Note that the experimental challenge in these systems is 
not limited to trapping and cooling down all the species of 
the relevant multiplet, but doing so while conserving the 
SU(N) symmetry. This goal is yet to be achieved in the 
BCS-BEC crossover regime. This is due to the fact that 
the external magnetic field used to control the interaction 
between the particles breaks the SU(N) symmetry. 

Experiments with N = 2 cold atoms are used as a play- 
ground for condensed matter physics Q . Similarly, high- 
spin experiments may shed light on other fields of physics. 
The ST/ (3) case is relevant for color-superconductivity in 
QCD, where the three quark colors generate an SU(3) 
symmetry. The SU(4) case is relevant for nuclear physics, 
where four types of Fermions exist (the proton and neu- 
tron with their spin up or down), and form an approx- 
imate SU (4) symmetry. It is well known that pairing 
plays an important role in nuclear physics Q. Therefore, 
the remarkable advance and flexibility of the cold atoms 
experiments gives the opportunity of making experiments 
relevant to other fields of physics which by themselves are 
hard to control and measure. 

The generalization of superfluidity theory to the case of 
more then two Fermion species arose soon after the BCS 
paper was publishe d 01 , presenting the case of two band 
superconductivity [10j . Since then, various techniques 
were used to deal with the SU (N) case, see for example 
Ref. Considerable theoretical attention was dedi- 



cated to the SU(S) case (lj - H^ and the St/ (4) case [l5j . 
Recently, guided by the experimental development, the 
SU(6) case was also studied [16(. 

A gapless superfluid, where a superfluid with gapped 
excitation spectrum coexists with a gapless branch of 
excitation, is a well known phenomena that naturally 
emerge in the case of particle symmetry breaking. For 
example, when there is a difference in the density or the 
mass of the interacting species, the result is a mismatch 
in their Fermi surfaces and therefore a gapless superfluid 
emerges 17|. Here, in contrast, we study a gapless super- 



fluid in the symmetric case, where the SU (N) symmetry 
is unbroken and completely holds. 

In this letter we focus on S-wave pairing in systems 
with odd number of species N = 2n + 1. Starting with 
the N = 1 case, i.e. one particle species, we note that 
pairing is forbidden and only a gapless excitation can 
exist. For the N — 3 case, it was proven that a gapless 
excitation coexists with gapped superfluid excitation (l3| . 
For the general case it was claimed but not proven, that 
for any odd N = 2n + 1 the ground state of an SU(N) 
symmetric system is a gapless superfluid |14j . Here we 
prove this hypothesis to be true in general. 

Let us consider a gas of N species of Fermions, inter- 
acting via a contact S't/(7V)-symmetric pairing potential. 
This system can be described by the generalized Hubbard 
Hamiltonian. 
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where (c Q ) is the creation (annihilation) operator for 
an a type Fermion, /i is the chemical potential, and 
g > is the attractive contact pairing coefficient. Note 
that the contact interaction coupling constant g must be 
regularized. This is usually done through the relation 
m/47ra s = g^ 1 + X^^Cfc) -1 between g and the scatter- 
ing length a s , which is the relevant physical parameter. 

Using the Hubbard-Stratonovich transformation [l8j 
to decouple the interaction in the Cooper channel, one 



2 



gets the partition function. 

Z = J D(*,tf) J nfl(Ao^A^)e- s , (2) 
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The field operator \P = (ipi ip2 ■■■iPn) is a vector of coher- 
ent states, and A a p = — A^ a is the complex pairing field 
describing an S-wave Cooper pair consisting of Fermions 
of types a and /3. The total magnitude of the pairing 
fields is given by |A| 2 = J2 a <p |A a ^| 2 , and Q is the ex- 
tended Nambu-Gor'kov propagator. 

At this point let us assume that the pairing fields are 
fixed in time and space A a p (x, t) = A a p . This simplify- 
ing assumption allows us to write Q~ Y in a block diagonal 
form in the momentum-frequency (k, iu) space. Each of 
these blocks is a 2N x 2N matrix of the form 



l(tw-f fc ) D 



(4) 



where 1 is the N x N identity matrix, = k 2 /2m — /i, 
and D is the pairing matrix D a p = A a p. The corre- 
sponding 2N spinors are (k, ioS), k, — iw)), and the 
sum on W is limited to positive Matsubara frequencies 
only. 

Note that the vector (\& 'J') is not a standard Nambu 
spinor. In order to consider all possible N(N — 1) pairing 
types the dimension of Q^ 1 must be 2N x 2N, and there- 
fore the spinors are 2iV-vectors. This raises the problem 
of performing the integration over the N- vectors ^,^>, 
which in ([3]) seems not to be Gaussian anymore. How- 
ever in the block diagonal form ((4]) , the integration is still 
Gaussian, and hence can be done quite easily. 

To show that, let us assume that ( Ul J is an eigen- 



eigenvector 



vector of Q 1 with the corresponding eigenvalue iuj — 
Xi. One can easily verify that there exists a twin 

ioj + Xi. Hence, the eigendecomposition of Q~ 



with the corresponding eigenvalue 



given by Q^Q 1 Q = Mu) 
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, and A is a real diagonal ma- 
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trix composed of the eigenvalues Xi . Using this transfor- 
mation the integration over the Fermionic fields reads 



Y\_(iuj - Xi)(—ioj - Xi 



(5) 

where Z)(^ r , 1 J r ) stands for integration over the fields 
\l/(k, iu), *(k, iuj), *(-k, -iu), *(-k, -iu). Restoring 



the momenta and frequencies summation one finally gets 

n<,i W ,k(* w - x i)- 

Now that we have established the Grassmannian in- 
tegration, we shall prove that for odd N > 1, there is 
always a gapless branch of excitations, i.e. one of the 
eigenvalues of Q^ 1 is simply iuj — To show that, we 
recall Jacobi's theorem [19j that if M is a skew symmetric 
matrix of order 2n + 1 then the determinant of M van- 
ishes. Since is such a matrix, one can conclude that 
there exists a vector v for which D'v = 0, therefore, 



Consequently 



D 







and ( .?„ ) are twin eigenvectors of 



with eigenvalues iw — and ioi + respectively. 
These eigenvectors and eigenvalues correspond to a gap- 
less Fermi liquid excitation spectrum, thus proving our 
claim. 

As an example, we explicitly consider the SU (3) and 
SU(5) cases. 

In the case of three Fermion species, the gapless branch 
of excitations, has the corresponding eigenvector 
/ A 23 \ 

where v = — A13 and D^v = 0. The spectra 

v a 12 ; 

of the two gapped branches are Ek = ±y^f + |Ap. The 
mean field values of A a p can be found from the saddle 
point approximation dS/dA a p — yielding three identi- 
cal gap equations, 
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Note that the same gap equation is valid for the SU(2) 
case, resulting in the same magnitude of the (total) pair- 
ing gap and the same critical temperature. The density 
equation for the a Fermion species is, 
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where f a = |A^ 7 | 2 /|A| 2 . If the system is balanced, i.e. 
contains the same number of particles for each kind, we 
get |A 12 | = |A 23 | = |A 13 | and f a = 1/3. That is to say 
that one third of the system remains gapless while the 
residual two thirds are superfluid. Note that since the 
different A a p's are equal, none of the species is unpaired, 
but instead there is a collective gapless excitation. 

In the case of five Fermion species, the gapless branch 



of excitations has the corresponding eigenvector 
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where v 



V 



-m 
m 
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Va = A ( 3 7 A (5e - A^A^ + A^ e A 7<5 , 
and D'v = 0. The spectra of the four gapped branches 



are E k = ±<J$ + ±|A| 2 ± ±^/|A| 4 - 4|?7| 2 , where |?7| 2 = 

X^al 7 ?"! 2 - Note that the cases where |A| 4 = A\ri\ 2 and 
\rj\ 2 = will be of specific interest, as the former has two 
doubly degenerated gapped branches and the latter has 
three gapless branches and two gapped branches. 

To conclude, we have studied a system of iV-Fermion 
species with an attractive, on-site, SU(N) symmetric in- 
teraction. We have proved that for odd N, the ground 
state of the system is a gapless superfluid. Doing so, we 
have used a generalized Grassmannian Gaussian integra- 
tion. The key point in our proof is the fact that the 
determinant of the skew-symmetric pairing matrix van- 
ishes if its dimension is odd. As a result a gapless free 
Fermi mode appears. The cases N — 3 and N = 5 were 
solved analytically at the mean field level to demonstrate 
our general claim. We hope that future cold atoms ex- 
periments would confirm our predictions. 
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